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ABSTRACT 



A reinvestigation was conducted of an anomaly in the experimental 
angular correlation for the gamma-gamma cascade of Ni^O which was 
reported by Verser and Conway in 1960 and 1961, respectively. The measure- 
ment of this correlation was being used to test angular correlation 
equipment built at the US Naval Postgraduate School. Two new experimental 
determinations of the correlation function with the same equipment as used 
previously indicated good agreement with the theoretical function. 

Verser' s and Conway's experimentally-determined angular correlation 
functions were found to deviate less significantly from the theoretical 
function than they had reported since their standard deviations had failed 
to take into account the variations of the characteristics of the electronic 
equipment due solely to electronic causes. It is concluded that the anomaly 
reported by previous experimenters was statistical in nature. 

Programs for the CDC 1604 computer have been developed in FORTRAN 
language for complete data reduction accepting as input the total counts 
in each channel, the total coincident counts, and other experimental 
parameters. The program provides for the calculation of the coefficients 
for a cosine series and a Legendre polynomial series expansion of the 
least squares fit to the experimental points. Other programs were 
developed to calculate the corrections for non-point detectors to be 
applied to the theoretical coefficients, and to calculate coincidence 
resolving times from data obtained during a large number of runs. 

Major improvements were made in the equipment to be used in future 
angular correlation experiments, including the construction and testing 
of both a slow and a fast coincidence unit. 

The author wishes to express his appreciation to Associate Professor 
H. E. Handler who has guided the development of angular correlation 
measuring equipment at the US Naval Postgraduate School for the last 
three academic years. His recommendations, particularly in the area 
of data analysis, and his perusal of the draft of this paper resulted in 
a still more thorough reinvestigation of the results of the angular 
correlation studies reported by previous experimenters at this School. 

It is a pleasure to ackowledge the guidance and aid received from Lt Cmdr 
P. A. Phelps (USN) in comprehending the theoretical aspects of the basic 
experiment, in becoming familiar with the electronic equipment, and 
in writing and debugging the computer programs. Appreciation is expressed 
to Mr. M. A. Brillhart for his construction of the coincidence units, and 
for his assistance in developing the mixer circuitry. 
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INTRODUCTION 

A. Purpose of Angular Correlation Experiments 

In gamma-gamma angular correlation experiments, one observes 
the probability of emission of the second quantum of a cascade as 
a function of angle relative to the first quantum. In general, 
this probability is an anisotropic function of the angle. Such 
experiments, in conjunction with other nuclear studies, can lead 
to identification of the multipole orders of the emissions and 
the total angular momenta of the original, intermediate, and final 
states of the nuclide under investigation. ^1,2,3 The orienta- 
tion of the total angular momentum vector of the intermediate level 
is correlated with the direction of emission of each quantum. The 
anisotropy of the probability distribution pattern will be attenuated 
if the orientation of the total angular momentum vector of the inter- 
mediate level changes between the times of these two emissions. Thus, 
one must study a nuclide for which the intermediate energy level has 
a mean life which is short relative to the period of precession of 
the angular momentum vector of that level. Only if the total angu- 
lar momentum vector of the intermediate level does not appreciably 
change its orientation between the emissions of the two quanta can 
the probability distribution as a function of angle between the two 
emissions be interpreted meaningfully. 

Generally, no parity assignments can be made for the nuclear 
levels involved based solely on a directional angular correlation 
study such as described above, since the angular distribution of 
the emissions is the same for electric and magnetic multipole radi- 
ations of the same order. By means of a similar polarization- 
direction correlation experiment, wherein one measures the orientation 
of the polarization vector relative to the plane of the two 

quanta, one can determine whether the emissions are electric or mag- 
netic multipoles. [ 4 , 5 j This information, coupled with known mul- 

tipole orders obtained from an angular correlation study, enables 
the assignment of relative parities to the three nuclear energy 
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levels. Absolute parity assignments require additional data* 

In the case of Ni^, an even-even nuclide, an even parity 
and zero total angular momentum are predicted for the ground 
state on the basis of the nuclear shell model, and are corrobo- 
rated by other experimental techniques unrelated to this paper . 
Studies of the decay scheme have shown that the spins of the 
first and second excited levels are 2 and 4 respectively . Cor- 
relation studies are consistent with this and demonstrate that 

these radiations are pure electric quadrupoles. Thus the com- 

60 

plete decay scheme for Ni is 



Decay Scheme of Ni 
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Since this decay scheme is well understood, this gamma-gamma 
cascade of Ni^ is frequently used to determine whether angular 
correlation equipment is functioning properly*, f 6 j It is used 
in this manner in this paper. 

B. Methods of Conducting the Experiment 

Since nuclei are generally randomly oriented, one cannot 
measure the radiation pattern of a given transition. One way of 
minimizing the random orientation of the nuclei is to apply a 
strong magnetic field or a strong electric field gradient at very 
low temperatures; radiation patterns have been measured using this 
technique. In the event of a cascade, the effects of this random 
orientation can be partially circumvented by establishing some re- 
ference direction with respect to which only certain of the second 
of the two radiations are measured. In the experiment discussed in 
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this paper, nuclei have not been aligned; the reference direction 

for observation is the direction of emission of the first quantum 

60 

in the two-gamma cascade of Ni , Since the lifetime of the intern- 
mediate nuclear level for this nuclide is short in comparison with 
the period of precession of its angular momentum vector, the ori- 
entation of the angular momentum vector does not vary significantly 
between the times of emission of the two quanta • Consequently, an 
unattenuated angular correlation between the directions of emission 
of the two successive radiations is observed. 

C. Interpretation of the Experimental Data 

By measuring the total , coincidence counting rate at several 
angles 9. between successive quanta, and applying corrections for 
accidental coincidences and for the attenuation of the correlation 
function due to the use of detectors of finite size, 7,8* one can 
determine the relative likelihood of emission of the second quantum 
at these angles relative to the direction of emission of the first a 
The experimental estimate of the correlation function at a specific 
angle 0^ is calculated according to the formula 



W (6.) 
ex i' 



N 1 N 2 



( 1 ) 



where C is the genuine coincidence counting rate at the given angle, 
G 

and Nj and ^ are the individual channel counting rates at that 

angle corrected for the system dead time and background. The least 

squares fit to these data is the experimental correlation function 

Wi s (9). 7 W^ s (0) ex P resse d in terms of a linear combination 

of even powers of cos 9 or in terms of a linear combination of even 

order Legendre polynomials. In the case of the gamma gamma cascade 
60 

of Ni , the theoretical correlation function, W^^CO) 3 corrected for 

the solid angles subtended by the detectors used in this experiment 

* 

is given by 



w th (e> = 1.00000 + 0.09960 P 2 (cos 0)+ 0.00838 P 4 (cos 0) . (2) 
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Experimental results can be compared with theory by comparing the 
coefficients calculated in the least squares fit with the above 
coefficients. In addition, one can compare the experimentally- 
determined anisotropy with the theoretical value. The anisotropy 
is defined as follows. 



W ls (180°) 
“ W ls (90°) 



( 3 ) 



The theoretical value of R, corrected for the geometry of this 
experiment is + 0. 16221 . A close agreement between the experi- 
mental and theoretical coefficients and between the experimental 
and theoretical anisotropies provides a check of proper equipment 
operation* 

Do Previous Work in This Field at U.So Naval Postgraduate School 
The basic equipment for this experiment was assembled in 
academic year 1959-60 by Verser [ 9 , who conducted an angular 
correlation experiment in the 90-180° quadrant and reported agree- 
ment with the theoretical anisotropy but reported marked deviation 
from theory in the shape of the experimental curve. The equipment 
was modified slightly in academic year 1960-61 by Conway £loj who 
performed the same experiment in the same quadrant and also reported 
a similar deviation from theory in the shape of the experimental 
curve. Fig (2) shows the results of these investigations. Verser 
reported a dip in his experimental curve which deviated by more than 
two standard deviations from the theoretical curve in the 120-135° 
region. Conway reported a similar dip in the vicinity of 135°, and 
a less pronounced dip in the 120° region, 

Conway wrote two programs in CDC 1604 computer assembly routine 
language, one for the least squares fit of a linear combination of 
Legendre polynomials of even order to the experimental data, and a 
second to calculate the solid angle corrections necessitated by the 
finite size of the detectors. The least squares program required a 
considerable amount of hand-calculated data reduction before it could 
be used. 
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Possible mechanisms underlying the reported anomalous dip 
in the experimental curve as suggested by Verser and Conway in~ 
elude the following: 

1, A change in the conditicnof the source, possibly 
crystallization resulting in preferential align- 
ment of the nuclei within the crystal due to in- 
ternal electric fields . 

2, Fluctuations of ambient temperature which resulted 
in changes in the operating characteristics of the 
equipment „ 

3, Uncompensated variations in line voltage, 

4, Perturbing effects of external magnetic fields, 

5, Improper orientation of the source relative to the 
detectors . 

It was recognized that the likelihood of having any noticeable 
effects due to Items "l" or ,, 4" above at room temperature was ex- 
tremely small since the Boltzmann factor for the relative popula- 
tions of the energy levels is nearly unity except under the influ- 
ence of extremely strong fields. 

In addition, the possibility was recognized of there being in 
fact no real anomaly, but that the discrepancy was simply an un- 
likely statistical occurrence or that some systematic error existed 
in the experimental technique, in the processing of data by hand, or 
in one of the previous computer programs. The likelihood of this 
deviation from the theoretical being simply a statistical fluke be- 
came greater when it was recognized that the variation in the re- 
solving time of the coincidence unit was not taken into account in 
computation of the standard deviation of the correlation function 
evaluated at the angles at which the experiments were conducted, 

E, Scope of the Investigation 

Many experimenters have reported agreement with the theoretical 
gamma-gamma angular correlation function for Ni^ j^ll,12j, yet two 
experimenters at this School independently reported a significant 
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deviation from theory in the 105-135° region. Thus it was clear 
that before one could use this equipment for the experiments for 
which it was assembled; namely, experiments with short-lived nu- 
clides, further investigation was necessary either to arrive at 
a mechanism for the reported anomaly, or to demonstrate that the 
anomaly, under more closely controlled experimental conditions, 

is not observed. In the investigation of this anomaly, certain 

1 

improvements were made in the experimental equipment; new pro- 

i 

grams in FORTRAN language were developed to compute the solid 
angle corrections to the angular correlation function and cal- 
culate the least squares fit; and experiments were conducted 
under closely controlled conditions in both the 90-180° and the 
180-270° quadrants to see whether a similar dip would be observed. 

Two additional coincidence units were constructed. One of 
these is a slow coincidence unit which, unlike that used by Verser 
and Conway, has a resolving time which is independent of counting 
rate. This characteristic will be required in conjunction with 
future experiments with short-lived nuclides. The second coinci- 
dence unit has a resolving time an order of magnitude shorter than 
that of the slow units. The use of the fast unit in conjunction 
with a slow unit enables one to achieve better precision by reducing 
the accidental coincidence counting rate. 
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REINVESTIGATION OF THE EARLIER REPORTED ANOMALY 

A. Equipment Used 

The experimental setup employed by Verser and Conway is des- 
cribed in detail elsewhere |^9,loj. Figure (3) indicates the general 
test setup and Figure (4) shows the angular correlation table with 
the source, crystals, photomultiplier tubes, and preamplifiers. The 
same slow coincidence unit as used by Verser and Conway was used in 
the angular correlation experiments in the 90-180° and in the 180- 
270° quadrants. The only significant change in the setup for these 
experiments consisted of adding a voltage regulator to minimize the 
effects of uncompensated line voltage fluctuations . Regulated vol- 
tage was then supplied to those parts of the electronic apparatus 
most likely to distort the experimental measurements when under the 
influence of line voltage fluctuations. These included the coinci- 
dence unit, both linear amplifiers, and the high voltage power supply 
for the photomultiplier tubes. 

B. Coincidence Resolving Time Tests 

The accidental coincidence counting rate is calculated from the 
usual formula 



C A = 2 r NjN 2 , (4) 

where ^ is the coincidence resolving time. Before and during the 
angular correlation experiments, scores of tests were conducted to 
determine the stability and the magnitude of T . The unit exhibited 
random fluctuations in resolving time as much as three times as large 
as those that can be explained in terms of nuclear statistics. The 

r ClD 

fractional variation in ' , y was observed to be approximately 
twice that reported by Verser and Conway. The measured mean is 
0.2470 microseconds. The variance of a single measurement of T is 
0.0044 microseconds; this includes the variation due to non-nuclear 
fluctuations and was used in calculating the variance of the acci- 
dental coincidence counting rates in this work. Since the resolving- 
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Figure 3, Experimental Setup of Verser and Conway 




Xtal* s 
PM' s 

Amp 1 s 



Coinc 



= Harshaw Scintillator Crystals, Type 8S8/2,(2 n diam x2 n thick) 

= RCA 6342 Photomultipler Tubes 

Power for PM tubes from Hamner H.V\ Power Supply Model N-401 

- Hamner N302 Amplifiers - each includes a differential pulse 
height analyzer which supplies a standard pulse out ("PHS" out) 
(A linearly amplified input is also supplied as an output but 
was not used in this experimental setup) . 

= Coincidence Analyzer * Atomic Model 502A. 



Scalers - Baird-Atomic, Model 1283 
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Fig. 4 Angular Correlation Table 



time characteristics were investigated and documented far more com 
pletely than in previous years, the observed instability was felt 
to be similar to that which would have been found by the earlier 
experimenters had they investigated this problem more thoroughly. 

An inclusion of the effect of this fluctuation in the analysis of 
their data would decrease their reported precision. 

C. Daily Checkout Procedure and Instrument Log 

A daily checkout procedure was devised to ensure that exactly 
the same checks of the operating characteristics of the electronic 
equipment were performed each day in the same sequence and with the 
equipment in the same operating configuration at each stage of the 
checkout . 

v An instrument log was developed in which the experimental 
equipment layout was recorded, and monitoring points for key wave 
shapes were indicated. Photographs of these key wave shapes are 
included, as well as the record of the oscilloscope settings used 
when the photographs were taken; a sample page from this log is 
shown in Figure 5 ° This log has proven useful in trouble-shooting 
the extensive circuitry and for monitoring changes in operating 
characteristics, such as pulse rise times, which could degrade system 
performance . 

D. Scaler Dead Time 

As a part of the complete checkout and alignment of all elec- 
tronic gear prior to conducting the angular correlation experiments, 
two-source experiments were run to determine the system resolving 
time (primarily scaler dead time) . Several of these measurements all 
led to the same surprising result, namely, a negative dead time for 
the system. Further investigation of this problem involved supplying 
pairs of pulses from a double pulse generator at a constant repetition 
rate and varying the time between the pulses of each pair. Expected 
behavior was that with the 1000 per second repetition rate for the pair 
of pulses, 2000 counts per second would be registered on the scalers as 
long as the separation between the pulses of each pair exceeded the 
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Monitor Point 15 Monitor Point 22 
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Fig. 5. Instrument Log 





resolving time of the system. When the separation between the two 
pulses was less than the resolving time of the system, 1000 counts 
per second would be expected. What was observed was the recording 
of 2000 counts per second when the interval between the pulses of 
each pair was greater than about 5-10 microseconds, then as the in- 
terval was decreased, the count rate recorded by the scalers in - 
creased , on occasion by as much as 257 e ; i.e. , to about 2500 counts 
per second. 

Since the experimentally-determined correlation function is 
not sensitive to the introduction of these spurious counts because 
the error introduced into the single channel counting rates is es- 
sentially constant for all angles and the effect on the coincidence 
counting rate is negligible, this anomalous behavior had no signifi- 
cant influence on the experimental results reported herein. New 
scalers will be employed for future experiments. 



E. 



Calculation of the Standard Deviation of W (0.) 

ex i/ 

The method of calculating the standard deviation of the experi- 
mental estimate of the correlation function has been the subject of 
extensive consideration. The experimental estimate of the correlation 
function at a specific angle 9^ is calculated according to the formula 



W (0.) = 
ex N 



N 1 N 2 



( 1 ) 



where C is the genuine coincidence counting rate at the given angle, 
G 

and and are the individual channel counting rates at that 
angle corrected for the system dead time and background. The genuine 
coincidence counting rate is calculated from 




where C = the genuine coincidence counting rate, C_ - the total coinci- 
G o T 

dence counting rate, and = the accidental coincidence counting rate. 
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C T and the single channel counts were measured directly by scalers 
and the accidental coincidence rate was calculated from 



C A = 2 T N 1 N 2 • (6) 

It can be seen that the standard deviation of W (9.) depends 

ex' l 

almost entirely upon that of C . The standard deviation of C is 

G G 

obtained by taking the square root of the sum of the variances of 
the two presumably independent quantities, and 



(T 



(C G> “ \ 



<r 2 (c T ) +<t 2 (c a ) 



(7) 



Verser and Conway calculated this standard deviation assum- 
ing that the distributions of and are Poisson in nature and 
therefore that the variances are expressible as: 

(T 2 (C a )=^ and (T*<C t ) = (8a, b) 

where t is the duration of the measurement. Therefore they calcu- 
lated (C ) from the relation 
G 



r(6 G ) - 




(9) 



If the equipment used in the experiment were perfectly stable, with 

no changes in T due to changes in electronic characteristics, then 

all fluctuations in C and C would be consistent with Poisson dis- 

A G 

tributions and these equations would apply. As pointed out in 
2 ’ 

Section B, ff~ (C ) when calculated from equation 8a is too small. 

To calculate a realistic variance of C^, one must provide not only 
for the statistical variations due to the random nature of the radio- 
active decay process, but also for the variations in the resolving 
time of the coincidence unit due to changes in its electronic oper- 
ating characteristics. 
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If electronic instabilities were present and nuclear fluc- 
tuations were small in comparison thereto, then the variance of 
would be described adequately by the relation 









cr) 



-i£ 



r 



b 



lj e< ^ u iov\ ( £) 



Then 





( / o) 



where (T~ _ denotes the standard deviation due solely to the elec- 
elec 

tronic instability of the coincidence unit. 

Since under experimental condit ionsneither of the two limiting 

conditions obtains; i.e., perfect equipment stability in which case 

the variance of would be Poisson, or infinite counting time in 

which case nuclear fluctuations would be insignificant in comparison 

with (P (C ), both statistical deviations must be taken into 

elec A / 2 • 

account in calculating (C ) . While the technique for the rigorous 

A 

statistical combination of these variances is not certain in spite 
of efforts to arrive at a valid analytical expression, the rather 
naive approach of adding them would seem to be a conservative ap- 
proach which applies to the experimental situation. Thus the ex- 

2 • 

pression used in this report for calculating (p (C ) is: 

G 




It can be seen that this expression gives the correct standard 
deviation in the two limiting cases discussed earlier, and it is 
assumed that the error resulting from applying this expression in 
the actual situation is small. 
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On the basis of the experimental determination of T , the third 
term of this expression is of the same order as the first two; hence, 
failing to take it into account results in unrealistically small esti- 
mates of the uncertainties in the experiment. 

F 0 Preparation and Check of New Computer Programs 

Conway prepared two programs in 1960-61 in 1604 computer assem- 
bly routine language. The purpose of one was to calculate the solid 
angle correction following the scheme of Rose j^8 j to apply to the 
theoretical correlation function to account for the attenuation re- 
suiting from the finite dimensions of the detectors. The second 
program calculated the coefficients for a least squares fit £ 7 j of 

even order Legendre polynomials to W (0.). 

ex l 

As mentioned previously, difficulties were experienced in at- 
tempting to run these programs. In addition, the least squares pro- 
gram required the following prior manual computations 

1. application of scaler dead time corrections 

2. correction of counting rates for background 

2 • 

3. calculation of C A and (T (C 4 ) 

A A 

. 2 9 

4. calculation of C and (p (C ) 

G G 

r\ _ 

5. calculation of W (9.) and <p [w (Q.)| 

ex' V i ex 1 J 

6„ calculation of a weighted mean of ^ ex (®^) at each angle 

7. calculation of the standard deviation of the weighted 
mean at each angle 

The values of the weighted mean of W^CO^) and the standard deviation 
of the weighted mean at each angle were the entering arguments for 
Conway's least squares program. 

The original intent was to program the above listed preliminary 
computations only, and to append Conway's least squares program as a 
subroutine. It was not possible to append his program to the FORTRAN 
program, however, due to the incompatibility of the two computer 
languages. A FORTRAN program was written to accept as input the total 
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counts registered by the scalers and thus save considerable time 
and achieve greater accuracy by avoiding the errors likely to creep 
into the lengthy manual computations. 

Separate programs were written which would compute the coeffici- 
ents for the even power cosine series as well as the coefficients for 
the even order Legendre polynomial series representation of W^(9), 
calculate the solid angle corrections to be applied to W^(9) , and 
calculate T from the data obtained during coincidence resolving 
time runs. 

The above programs provided independent checks of Conway's 
least squares and solid angle correction programs. The new least 
squares program proved extremely useful during the reanalysis of 
earlier data as it permitted repeated reduction of the same experi- 
mental data with incrementally varied input parameters. 

The FORTRAN programs together with a glossary of terms is in- 
cluded in the Appendix fA thru E.^ 

To check the operation of the FORTRAN programs, they were used 
to perform the same least squares and solid angle correction computa- 
tions as performed earlier by Verser and Conway, and the results com- 
pared with their earlier results, as follows. 

1. The same geometrical parameters as used by Verser and 
Conway were used as the entering arguments for the 
FORTRAN solid angle correction program. Good agree- 
ment was reached with Verser* s hand-calculated values 
and Conway's computed values. 

2. The raw data (scaler counts) or Verser were used as 

the entering arguments for the FORTRAN least squares 

program. Good agreement was reached with Verser 1 s 

values of the mean values of W (9.), and with his 

ex x 5 

coefficients for the Legendre series representation 
of W^ g (9) thus indicating proper operation of the 
preliminary and subsequent portions of this program. 
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3. The results of analysis of Conway* s data with the 
FORTRAN program did not agree with the results of 
analysis which he reported. Some significant 
arithmetical errors were uncovered in the hand 
calculations which he was required to perform 
prior to entering his program. When these errors 
were corrected, agreement was realized in the re- 
sults of analysis by the two programs, providing 
another check of the FORTRAN least squares program. 

G. Reanalysis of Earlier Data 

Once assured that the new FORTRAN programs were operating pro- 
perly s Verser's data and a portion of Conway 1 s data were reanalyzed, 
including in the calculation of the variances the effect of the 
electronic instability of the coincidence unit. 

In checking over the data collected by Conway, it was noted that 
they were collected during two relatively short periods of time spaced 

nearly four months apart. During these four months a significant change 

60 

in the strength of the Co source (T ~5.2 years) occurred; this 
had not been corrected for in his data analysis. The reanalysis of 
his experiment was done using only his later set of data. 

It became apparent that the scatter of Conway's data about W^ g (9) 
was greater than the standard deviations would imply. Consequently 9 
estimates of the values of (T^W ex (0^)j were calculated based upon the 
actual residuals between the experimental points and Wj g (9) for his 
later set of data. This calculation, which took into account the elec- 
tronic variations of the equipment, resulted in values for 




W 

ex 

W 



ex 




(V 



approximately three times as large as the values which Conway quoted. 

The results of the reanalysis of the previous data are tabulated 
below and are presented graphically in Figure 6a and 6b. It is evident 
that the anomaly in the 120-135° region still exists in Verser's results, 
though it is not as significant statistically as previously reported due 
to the larger standard deviations. The analysis of Conway's data now 
shows no evidence of a dip. 
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Reanalysis of Data of Verser and Conway 
Fig. 6a, b 
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Table 1 

Results of Reanalysis of Earlier Data 





Legen 

a 

o 


dre Coefficients 
a 2 


a 4 


Theoretical 


1.0000 


0.0996 


0.0084 


Verser 


1.0000 + 0.0042 


0.0960 + 0.0068 


0.0337 + 0.0087 


Conway 


1.000 + 0.011 


0.126 +0.020 


0.032 + 0.023 



An: 


Lsotropy 


Theoretical 


0.162 


Verser 


0.171 ± 0.026 


Conway 


0.220 + 0.079 



The possibility was investigated that errors could have been intro- 
duced into Verser 1 s data analysis through the use of incorrect parameters. 
In investigating this possibility, raw data gathered by both Verser and 
Conway were used as the entering arguments for the FORTRAN least squares 
program, and these data were reduced while systematically varying the ex- 
perimental parameters to see what effect this would have on The 

results of this analysis are as follows: 



Table 2 

Effect on W^CG) of Varying Input Parameters 



Parameter Varied 


7o Variation in 


Resulting % Variation 




Parameter 


Legendre Coefficients 


Anisotropy 






a 2 


a 4 




r 


+ 10% 


+8% 


+8% 


+8% 


Background 


+ 6% 


0% 


0% 


0% 


Scaler Resolving Time 


+500% 


+3% 


+3% 


+5% 



It can be seen from the above that the shape of W^ s (Q) is relatively in- 
sensitive to changes in background and changes in the resolving time 
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of the scalers. However the shape of W.. (0) varies appreciably 

l s 

when T varies. The reason for these changes in W^ s (9) is clear 

from the following. If the numerical value used for T is greater 

than the actual T of the coincidence unit during the experiment, 

then too many accidental coincidences are subtracted from the total 

coincidences in calculating the number of genuine coincidences. 

Thus the apparent value of C is smaller than was actually realized 

G 

experimentally. Since the error introduced into is constant at 

all angles, assuming that and are independent of 9, which is 

approximately true, the error introduced into C is also constant 

G 

at all angles. Since in these experiments, the ratio of the acci- 
dental coincidence counting rate to the genuine coincidence counting 
rate is about 2:3 and it can be seen from the expression 

W (180) \ 

rar - 1 - • 167 (3) 

Is 




that W- (90) is more sensitive than W. (180) to changes in the genuine 
Is Is 

coincidence counting rate, it is apparent that the calculated anisotropy 
increases when the coincidence resolving time is increased. It is 
stressed that the shape of W^ g (0) is distorted by using a value of T 
in the data reduction which does not correspond to the value of T 
actually characteristic of the coincidence unit during the experi- 
ment. The shape of W.. (9) is independent of T as long as the true 

Is 0 

experimental value of T is used in calculating C . 

A 

The effect upon the shape of W.. (9) of using incorrect para- 

x s 

meters for the geometry of the experiment was investigated also, by 
progressively varying the entering arguments for the FORTRAN solid 
angle correction program. The results of this analysis are given in 
in Table 3. 
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Table 3 

Effect on Solid Angle Corrections of Varying Input Parameters 



Parameter Varied 


Amount Varied 


Relative Ch 
Absorption 


ange in 
Coefficients 


Linear absorption 
coefficient of the 
crystal 


+ 107, 


0.00% 


0 . 00% 


Distance from source 
to crystal 


+ 8.5% 


0.00% 


0 . 00% 


Crystal radius 


+ 1.0% 


0 . 00% 


0 00% 


Crystal thickness 


+ 1.0% 


0 . 00% 


0.00% 



It can be seen that the shape of W^ s (9) is relatively insensitive to 
changes in these geometrical parameters over the range investigated. 

It is concluded from these investigations of the effects upon the 
least squares analysis and solid angle correction computations of varying 
the entering arguments, that none of these can account for the dip re- 
ported by Verser, 
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H Analysis of New Data 



The results of the analysis of the data obtained with the old 
slow coincidence unit in experiments in two separate quadrants are 
indicated graphically in Fig, 7. When the data from each of these 
two experiments were reduced separately, W^ g (9) a g re ^d with W t ^(9) 
within the statistical error of the experiment. Since the angular 
correlation function is symmetrical about 180°, it was possible to 
combine the results of the analysis of the data obtained during the 
experiment in the 90 - 180° quadrant with those obtained during the 
the experiment in the 180 - 270° quadrant. The purpose of combining 
the results in this way was to benefit from the decrease in (f W ex (£K) 
resulting from the analysis of a greater quantity of data, thus being 
able to compare more precisely W^(9) a nd W t ^(9) 

Due to the time interval between the two experiments, during 

which the source decayed significantly, the raw data from the two 

experiments cannot be combined. Consequently, the data obtained from 

each experiment were analyzed independently, then the normalized 

values of the weighted mean of ^ x (9^) from the two experiments were 

combined using as weighting factors the reciprocals of their variances. 

These composite values of W (9.) were then used as the entering 

ex i 

arguments for the least squares subroutine to generate a new (9) . 

Is 

The results of the new experiments are tabulated below. 



Table 4 

Results of Analysis of New Data 





i 

Experimental 


Theoretical 


Quadrant 


Composite 


180°-270° 


90-180° 


Anisotropy 


0.151 + 0.029 


0.165 + 0.034 


0.158 + 0.018 


0.162 


Legendre 










Coefficients 










a 

o 


1.0000 + 0.0057 


1.0000 + 0.0057 


1.0000 + 0.0040 


1.0000 


a 2 


0.092 + 0.010 


0.1062 + 0.0094 


0.0987 + 0.0069 


0.0996 


a. 


0.010 + 0.012 


-0.005 + 0.012 


0.0036 + 0.0085 


0.0084 


4 


— 






i 
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Results of Analysis of New Data 
Fig. 7. 
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I„ Comparison of Experimental Results with Theory 

All of Verser*s data, and the portion of Conway* s data used in 
the reanalysis of his data (see Section G) were combined with the 
data obtained during the two new experiments according to the scheme 
discussed in the previous section. That is, the data from each of 
the four experiments were analyzed independently, then the weighted 
means of the normalized mean values of W ex (9..) from the four experi- 
ments were used as the entering arguments for the least squares sub- 
routine. The results of this combination of these data are indicated 
graphically in fig* 8, and are tabulated in Table 5» No significant 
deviation from W t ^(9) is evident. 

While the real demonstration of the proper functioning of the 
experimental equipment lies in the results of the new experiments in 
two quadrants, it is of interest that no significant deviation from 
W t ^(9) occurs even when Verser's and Conway* s results are combined 
with these. It is concluded that the dip which they reported is 
strictly a statistical fluke and that there is no actual physical 
effect involved. 



Table 5 



Results of Composite Data Analysis 
(Composite Results based upon 4 experiments) 



, 


Experimental 


f 

Theoretical 

_ _ ) 


Anisotropy 


t- • ~~ • 

0.169 + 0.018 


0.162 


Legendre 






Coefficients 






a 


1.0000 + 0.0028 


1.0000 


o 






a^ 


0.1000 + 0.0047 


0.0996 


2 


— 




a. 


0.0183 + 0.0058 


0.0084 


4 
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Results of Composite Data Analysis 



Fig. 8. 
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DEVELOPMENT OF NEW INSTRUMENTATION 

A. New Slow Coincidence Unit 

While investigating the earlier reported anomaly in W^(Q) in 
the 105-135° region, work proceeded on building and testing a new 
slow coincidence unit incorporating certain portions of the design 
reported by Garg 13^ . This design features a variable resolving 
time which is independent of the counting rate. The slow coinci- 
dence unit used by Verser and Conway lacked both of these features 0 
The coincidence unit which was constructed differed markedly 
from Garg's only in that the separate channel amplitude discrimi- 
nators of Garg's unit were deleted since differential discrimina- 
tion was performed earlier (Fig. 9 ). Figure 10 gives the com- 
plete circuit diagram for the unit which was constructed. 

The new slow unit was tested by incorporating it into the 
circuitry used for the experiment in the 90-180° quadrant in 
parallel with the slow coincidence unit which Verser and Conway had 
used. The correlation function measured by the two slow coinci- 
dence units agreed within the uncertainties in their resolving times. 
The results of numerous resolving time tests of the new slow coinci- 
dence unit indicated a greater long time drift than is characteristic 
of the old slow unit. This increased drift of T necessitates more 
frequent determinations of its magnitude and required using different 
values for T during the angular correlation runs, rather than using 
a mean value of f based upon a large number of resolving time runs 
such as was done with the old slow coincidence unite 

B. Requirement for a Fast Coincidence Unit 
I. General 

A coincidence unit of significantly shorter resolving time 

than that of the unit used by Verser and Conway was required to 

improve the statistical error in the experimental estimate of 

the correlation function, W (9.). Better statistics were re- 

ex' i' 

quired for the reinvestigation of the earlier reported anomaly 
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Figure 9o Block Diagram, New Slow Coincidence Unit 
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Pulse Shapers Coincidence Amplifiers 



to figure 10: 
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Figure. 10 a. Circuit Diagram of New Slow Coincidence Unit 
Initial Stages 
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2 . 



in the value of W (0.) for Ni^ in the 105-135° region, and 

for future experiments with short-lived nuclides produced in 

the US Naval Postgraduate School Van de Graaff generator. 

It was not possible to reduce the resolving time of the 

slow coincidence unit below a value T . which is determined 

mm 

by the width of the window of the differential pulse height 
analyzer which precedes it. This time corresponds to the dif- 
ference in time between the crossing of the lower voltage dis- 
crimination level of a pulse of the lowest permissable amplitude 
in one channel and the time of crossing of the lower discrimi- 
nator level of a simultaneous pulse of the highest permissable 
voltage in the second channel. Because of this irreducible 
minimum value for T y and since a resolving time larger than 
7\ n gives too many accidental coincidences for an acceptable 
variance of the experimentally-determined correlation function, 
a different mode of operation was required. The scheme used 
consisted of adding a fast coincidence unit in parallel with 
the slow unit and mixing the output pulses from the two units; 
this scheme reduced the accidental coincidence counting rate 
because of the decreased T of the fast unit, yet retained the 
amplitude* discriminating characteristics of the differential 
pulse height analyzer employed with the slow unit (Figure 11 <>) 
The fast coincidence unit was constructed after a design 
by Draper and Fleischer (14) and was adjusted for a resolving 
time of approximately 20 nanoseconds, thus decreasing C by 

A 

about a factor of 10 from that obtained with the slow units. 
Relation between^ and 



It is desired to minimize 

r I 






( 1 ) 



W (0.) 
ex v i 7 



Since W (0 . ) = 
ex v i 



N 1 N 2 



and the fractional deviations of N, 
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Figure 11, Block Diagram of Experimental Setup 
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Slow Coincidences Atomic Unit 

Slow Coincidences New Unit 

Fast-slow Coincidences (fast 6c New Slow) 





and N 



V 



are very small compared to the fractional deviation of 



<T 



“ex< 6 l>1 . 



(12) 



w (0.) 

ex 1 / 



Effect on (T (C ) of Decreasing the Coincidence Resolving Time 
G 

The variance of C is calculated from equation H , 

G 



2 * T 

(T (C G ) = -f 



<n. 



ele c 

r 



(Jh 



• 2 
C A + 



t 



It is evident from the relation C = 2TN-,N~ (equation 6 1 that 

A 12| I 

decreasing the resolving time of the coincidence unit 7] by an 
order of magnitude results in a similar decrease in C^. This 
decrease is reflected as a decrease in the third term of equation 
11 by an order of magnitude since t is constant in these experi- 
ments o Assuming no change in the fractional resolving time fluctu 



mining 
,1 erf 



ation, | -jr- j 
by a factor of about 100 



the second term of equation 11 would decrease 
The genuine coincidence counting rate 



is independent of the coincidence unit resolving time, as long 
as the latter exceeds the average interval between the emissions 
of the two quanta; this is the case for all three coincidence 



units discussed in this report c 



Thus a decrease in C also re~ 

A 

This decrease in 



suits in a decrease in C since C = C + C , 

T £ T A G 

C results in a decrease in _T. Therefore all three terms of 

T t 

equation 11 decrease when T decreases. 

Relative Effect of Coincidence Resolving Time Drift 

2 • 

The previous section indicated how (j~(C ) is decreased by 

G 

decreasing T . Then from the relation 



<r 



t H e*<y] 



<T< c r> 



equation 12 , 
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it is evident, 



since is 



constant, 



that 



<r 



Kx<°l>] 

U C x<°l> 



has been decreased appreciably. 



The assumption that 



\<T i 

elec 

\ r 




did not vary was made in the last section to simplify the dis- 
cussion. It is evident, however, that due to the large decreases 

in the values of all other terms in equation 11 , including the 
* 2 

term which multiplies 
latter can be tolerated while still decreasing very significantly 
the fractional variance of W ex (CL) * Thus the electronic stability 
of the fast coincidence unit is not as critical as it would be 
with a unit of larger coincidence resolving time. 



^lec ^ ^ ) 

r 



, larger values of the 



Co Operation of the Fast Coincidence Unit 

The basic theory of operation of the fast unit and the circuitry 
for the Draper Fleischer unit are given in detail elsewhere ^ 14] and 
will not be repeated heref Either single channel becomes operative 
when the input pulse reaches 68% of its maximum amplitude, regardless 
of the value of that amplitude within the range 20 100 volts. The cir- 
cuit is so designed that resolving times on the order of 100 times 
shorter than the rise time of the input pulses can be achieved. This 
characteristic is required since the rise time of the pulses from the 
sodium iodide scintillation crystals is on the order of 200-300 nano- 
seconds o 

Only the coincidences of 1.17 MEV and 1„33 MEV gamma rays are of 
interest in this experiment with Ni^. Amplitude discrimination which 
was performed by the differential pulse height analyzers which preceded 
the slow coincidence unit selected the pulses corresponding to the 
photopeaks of the two gamma rays. All coincidences, regardless of the 
pulse height, were detected by the fast unit since no amplitude dis- 
crimination preceded it. The output pulses of the slow coincidence 
unit and those of the fast unit were passed through a mixer circuit 
which effectively functioned as another coincidence unite An output 



* A block diagram of the fast unit is provided in Fig. 12. 
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Block Diagram, Fast Coincidence Unit 
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After Bell 6c Petch 6c 



pulse was obtained from the mixer only if there was, in a time interval 
of about 35 microseconds, a 20 nanosecond coincidence between pulses 
of any height a$ ascertained by the fast unit and a 200-nanosecond 
coincidence between pulses corresponding to the photopeaks of interest 
as determined by the slow unit. Thus by combining these two units 
in this manner, an energy-selective counting system with a 20-nano- 
second coincidence resolving time which was independent of the counting 
rate was obtained. 

Certain modifications were necessary to the fast unit and are 
given in Figure 13, The last two stages of the Draper-F leischer 
unit could not be made to function properly. This portion of the 
circuitry was therefore taken out, and a cathode follower was sub- 
stituted . 

The mixer circuit which contained circuitry for amplitude dis- 
criminating, shaping, amplifying, and mixing the fast coincidence 
pulse with the slow coincidence pulse was built up and tested stage- 
by-stage. The discriminator of the mixer unit is a uni-stable multi- 
vibrator with Zener diode stabilization. The inverter stage in the 
fast coincidence channel of the mixer results in a negative pulse 
of approximately 35-microsecond duration, the duration varying with 
the discriminator setting. This pulse is then mixed with a 0.7 micro- 
second negative pulse from the slow coincidence unit. Discrimination 
of the mixed output is performed by the scaler counting these fast- 
slow coincidences. The design of the mixer circuit is given in Fig. 14. 

Do Testing and Alignment of the Fast Coincidence Unit 

The alignment of the fast coincidence unit consisted of two 
major portions-achieving linearity in each channel, and assuring the 
detection of true coincidences by the coincidence stage. 

Draper and Fleischer ^lJ describe a means of testing for 
linearity; i.e., ensuring that the fast discriminator fires when the 
input reaches the same predetermined fraction of its full amplitude, 
over a wide range of amplitudes. They prescribe a mercury-relay pulser 
delivering a five-volt negative step with a rise time of seven nano- 
seconds which is split, one portion triggering a 10-mc ringing cir- 
cuit which provides the time reference, the other going through an 
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Discriminating Circuitry 
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Figure 14. Circuit Diagram of Mixer Circuit 



attenuator, linear pulse amplifier, and the coincidence unit 
(Figure 15 -a) . 

Several modifications of this test scheme were tried, due to 
the nonavailability of a pulse generator which could deliver a pulse 
of sufficiently large amplitude and short rise time to trigger the 
ringing circuit satisfactorily, and due to the lack of a compensated 
attenuator. These modifications met with varying degrees of success. 
When a double pulse generator was used as a signal source, and the 
amplitudes of the inputs to the separate channels of the fast unit 
were varied by adjusting the gains of the linear amplifiers, dis- 
tortion of the wave shapes resulted when the coarse gain settings 
were varied. This made it impractical to vary the amplitudes of 
the coincidence unit inputs by varying the gains of the linear 
amplifiers since the time reference for determining the delay within 
the coincidence unit itself was thus lost. Still greater wave shape 
distortion resulted from varying the gain of the double pulser. 

The procedure which has finally emerged as workable is outlined 
in Figure 15-b and consists of triggering the oscilloscope with 
the output of the fast unit discriminators and viewing the relative 
delay of nuclear pulses of varied amplitudes. The problem of wave 
shape distortion when the coarse gain settings of the linear ampli- 
fiers were adjusted has been avoided by limiting the adjustment of 
delay-equivalence to the 10-volt range of amplitudes used experi- 
mentally. Since a 10-volt window is used in the differential pulse 
height analyzer experimentally, an equivalent time delay over this 
range is all that is required. The inputs to the fast unit are varied 
over this 10-volt range by adjusting only the fine gain controls, 
avoiding the distortion of wave shape which accompanies changes of the 
coarse gain settings. The adjustment of the fine gain controls can 
be avoided by providing a nuclear source or a combination of sources 
with different decay energies which results in different amplitude 
signals being received by the fast unit simultaneously. The relative 
delay introduced into these pulses of different amplitudes is then 
observable directly. 
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Fig. 15. Fast Coincidence Unit Test Scheme 
(As described by Draper and Fleischer) 
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To scope input 



An alternate means of accomplishing this alignment consists of 
using test pulses from a double pulser, triggering the oscilloscope 
with the trigger output pulse of the double pulser and viewing the 
output of the fast discriminator of the fast coincidence unit in 
each channel. Again, the input amplitude is varied by changing the 
settings of the fine gain controls on the linear amplifiers. 

The procedures already described ensure only that the delay 
which a pulse experiences as it passes through the time compensating 
circuit of one channel is independent of amplitude „ The other major 
alignment requirement, which cannot be met until each channel has 
been aligned separately, is the assurance that the true coincidences 
are being detected by the unit. This requires identical total delay 
from the time of emission from the source through all the circuitry 
associated with each of the two channels and through the fast coinci- 
dence unit itself up to the coincidence stage. 

To determine the ability of the fast unit to detect true coinci- 
dences before it was completely aligned, it was used in parallel 
with a slow unit in the angular correlation experiment in the 90-180° 
quadrant. While it generally tracked the slow unit, it did so errati- 
cally. 

Variable delay lines have been installed to facilitate the align- 
ment for equivalent delay in the two signal channels to the coincidence 

stage of the fast unito The alignment is accomplished using the experi- 
60 

mental Co source with the experimental geometry and normalizing the 
total coincidence counting rate measured by the fast-slow coincidence 
detecting scheme to the genuine coincidence counting rate determined 
by the slow coincidence unit. After normalizing the total coincidence 
counting rate of the fast-slow combination in this manner, and making 
a graph of this versus the setting of the variable delay line in one 
of the time compensation channels of the fast unit, a plateau is ob 
tained the width of which corresponds to the resolving time of the 
fast coincidence unit. The delay line is then adjusted to the center 
of this plateau as determined from the graph. 
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Thus procedures have been developed for aligning the time com 
pensation circuit of each channel for linearity by means of either 
a nuclear source or a double pulser, and for aligning the two chan- 
nels of the fast unit for equivalent delay by means of a nuclear 
source . 

E Equipment Status 

The experimental equipment described herein including the new 
fast' slow coincidence detecting scheme is functioning well enough 
to be used in the experiments for which it was assembled; namely, 
the angular correlation of short' lived nuclides activated in the 
Van de Graaff generator. Due to limited unavoidable drifts in the 
electronic equipment, frequent checks must be made of operating 
characteristics, and appropriate adjustments made as necessary to 
ensure the closest possible control over the experimental environ- 
ment. The checkout procedures which have proved extremely useful 
in ensuring that easily' overlooked actions are not overlooked and 
that checks are made in the same sequence for the greatest possible 
reproducibility of procedures, should be continued 

Figure 16 illustrates the layout of the major electronic portion 
of the experimental equipment, all rack-mounted and mobile to facili- 
tate its use with or without the Van de Graaff generator. 
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Fig. 16 Electronic System 
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PROGRAM ANGCORR 

DIMENSION AN 1( 20 » 20 ) » AN2(20,20), CT(20,20), 83(20), 

K(20), AND 1 P ( 20, 20 ) , AND2 P ( 20 , 20 ) , CTD(20,20), CAD( 20,20) 
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CGD ( 1 1 , J J ) = C TO I I I , J J ) - C AD ( II » J J ) 

W(II,JJ) = ( CTO( I I , JJ )-CAD( I I , JJ) )/ ( (AND 1P(I I , JJ )-BACK1 ( 1 1 , JJ ) ) 
*(AND2P( II, JJ)-BACK2(II , JJ) )) 

SW( I I , JJ ) = Will , JJ)*SQRTF ( ( CTO ( II,JJ)/T(II,JJ) + (CAD(II,JJ)/T(II,J 
J ) ) + ( STAUCI I I , JJ >»CAD ( II , JJ ) / TAUC ( II,JJ))**2))/(CTD(II,JJ )-CAD 
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71 SUM3 = SUM 3+ (W( II,JJ)-B1 CJJ) ) * *2 
IF ( II-KI JJ) ) 72, 73, 73 
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1 CAO ( I 1 1 J ) 1 CGD ( I I , J ) , W(II,J), SWIII, J). JRUN ( 1 1 » J ) 
II = II+l 

IF III-K(JI) 85, 85, 87 
87 WRITE OUTPUT TAPE 4, 532 
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220 SUM5 = P0LY(K, I )»OMEGA(K»J )+SUM5 
225 PART l(I.J) = SUM5 
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ARG4F ( A » E f C ) = POL Y4F { C ) * S I NF ( C ) *E PSF ( A , E , C ) 
ARGP4F(A,B,C,D) = POL Y4F ( C ) *S I NF { C ) * EPS PF { A , B f C , 0 ) 



210 FORMAT (4F10.6) 

215 READ INPUT TAPE 3, 210, R ( 1 ) * TU), H(l), TAU(l) 
READ INPUT TAPE 3, 210, R(2), T{2), H(2), TAU ( 2 ) 
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260 SUMO = SUMO + ARGOF ( TAUI I ) , T(I), BP ( I ) ) 

SUMPO = ARGPOF ( TAU( I ) , R C I ) , BP(I), HU)) +4 . 0* ARGPO F { T AU ( I ) 
1RII), BP ( I ) +DELTAP { I ) , H(I)) ' 

BETA = BP ( I ) + 2.0*DELTAP( I ) 



270 SUMPO = SUMPO +2. 0*ARGP0F { TAU(I) , R(I), BETA, H(I)) +4.0*ARGPOF 
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Appendix E FORTRAN Program 
Glossary of Terms for the 
Least Squares Analysis 

1. Check of input parameters 

TAUC Coincidence unit resolving time (sec) 

STAUC Standard deviation of T (sec) 

TAUDl System dead time, Channel 1 (sec) 

TAUD2 System dead time, Channel 2 (sec) 

BACK1 Background, channel 1 (counts per sec) 

BACK2 Background, channel 2 (counts per sec) 

T Duration of run (sec) 

2 0 Check data input of N^N^and CT and input parameters 

Total counts, channel 1 

Total counts, channel 2 

CT Total coincidence counts 



3. Results of analysis 

ND0T1PRIME Count rate, channel 1, corrected for system dead 
time and background (counts per sec) 

Corresponds to ft., in text 
ND0T2PRIME Same, channel 2 1 



CTOTALDOT C T 

CACCDOT C A 
A 



CGENDOT 

W 



SIGMAW 



CL 



W (0.) - 

ex' i 



* 1*2 



/ i 



W (G.)l = 

ex i' 




X 


W (0 0 
ex i y 


c 


L J 


G 




- 



WBAR 



w_ (6i)- 

e x 








I 

^ 1 Iks* 




a 




where N is the number of runs at the angle 0^. 

EXPERIMENTAL FUNCTION W. (9) normalized to W n (90) and 

Is Is 

evaluated at 5° intervals. 

THEORETICAL FUNCTION W th (Q) normalized to W (90) and 

evaluated at 5° intervals. 



NORMALIZED EXPERIMENTAL WBARS AND SIGMAS 

WBAR W (0.) normalized to W (90) 
ex i J. s 



ANISOTROPY - P 



W. (180) 



1 



based upon Legendre series 
representation of 



ANISOTROPY - C 



Same 
of W 



, based on Cosine series representation 
Is ’ 



SIGMA Standard deviation of ANISOTROPY - P 
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